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ABSTRACT 

In the Cardassian model, dark energy density arises from modifications to the 
Friedmann equation, which becomes = g{pM), where g{pM) is a new function 
of the energy density. The universe is fiat, matter dominated, and accelerating. 
The distance redshift relation predictions of generalized Cardassian models can 
be very different from generic quintessence models, and can be differentiated with 
data from upcoming pencil beam surveys of Type la Supernovae such as SNAP. 
We have found the interesting result that, once VLm is known to 10% accuracy, 
SNAP will be able to determine the sign of the time dependence of the dark energy 
density. Knowledge of this sign (which is related to the weak energy condition) 
will provide a first discrimination between various cosmological models that fit 
the current observational data (cosmological constant, quintessence, Cardassian 
expansion). Further, we have performed Monte Carlo simulations to illustrate 
how well one can reproduce the form of the dark energy density with SNAP. 

To be concrete we study a class of two parameter (n,g) generalized Cardassian 
models that includes the original Cardassian model (parametrized by n only) as 
a special case. Examples are given of MP Cardassian models that fit current 
supernovae and CMB data, and prospects for differentiating between MP Car- 
dassian and other models in future data are discussed. We also note that some 
Cardassian models can satisfy the weak energy condition w > —1 even with a 
dark energy component that has an effective equation of state wx < — 1. 
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1. 



Introduction 



Recent observations of Type la Supernovae (Riess et al. 1998; Perlmutter et al. 1999) 
as well as concordance with other observations (including the microwave background and 
galaxy power spectra) indicate that the universe is accelerating. Many authors have explored 
a cosmological constant, a decaying vacuum energy (Freese et al. 1987; Peebles & Ratra 1988; 
Prieman et al 1995), quintessence (Wang & Steinhardt 1998; Caldwell, Dave, & Steinhardt 
1998; Huey et al. 1999), and gravitational leakage into extra dimensions (Deffayet 2001) as 
possible explanations for such an acceleration. Recently Freese & Lewis (2002) proposed 
Cardassian expansion as an explanation for acceleration which invokes no vacuum energy 
whatsoever.^ In this model the universe is flat and accelerating, and yet consists only of 
matter and radiation. 

In Cardassian models, the Priedmann equation is modified from H'^ — Snp/ (Sm^;) to 



where g{pM) is a different function of the energy density, pm contains only matter and 
radiation (no vacuum), H — a/ a is the Hubble constant (as a function of time), and a is the 
scale factor of the universe. Models with gravitational leakage into extra dimensions also 
give a modified Priedmann equation that can be cast into the form of Eq. (1), but with a 
vacuum component. 

The function g{pM) returns to the usual STrpM/i^fripi) during the early history of the 
universe, but takes a different form that drives an accelerated expansion after a redshift 
z ~ 1. Such modifications to the Priedmann equation may arise, e.g., as a consequence 
of our observable universe living as a 3-dimensional brane in a higher dimensional universe 
(Chung & Freese 2000). Alternatively, such a Priedmann equation may arise if there is dark 
matter with self-interactions characterized by negative pressure (Gondolo & Freese 2002). 

We wish to study the detectability of the altered Priedmann equations by upcoming 
observations of Type la Supernova such as SNAP. As discussed in Freese (2002), the redshift 
distance relationship predictions for generalized Cardassian models can be quite different 
from generalized quintessence models. It is the goal of this paper to see how well we can 
reproduce the correct form of this dark energy density in upcoming experiments. 

For concreteness, we investigate a particular version of generalized Cardassian cosmol- 
ogy. However, our results are intended to generalize to any Cardassian cosmology, i.e., to 



^The name Cardassian refers to a humanoid race in Star Trek whose goal is to accelerate expansion of 
their evil empire. This race looks alien to us and yet is made entirely of matter. 
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any function g{pM)- The particular model we study is meant to illustrate that, generically, 
modified Priedmann equations can lead to specific detectable predictions in experiments like 
SNAP. 



In this paper we study the following generalization of the original Cardassian model: 
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We call this model "Modified Polytropic Cardassian" (MP Cardassian)^ (Gondolo & Preese 
2002), where G — l/rripi is Newton's universal gravitation constant, pcard is a characteristic 
constant energy density, and where we take n < 2/3 and q > 0. The original power law 
Cardassian model corresponds to g = 1. 

Por comparison, we remind the reader of the original power law Cardassian model that 
was proposed in Preese & Lewis (2002), which had the following specific form of g{pM)'- 
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The first term inside the bracket of Eq.(2) and Eq.(4) dominates initially, so that or- 
dinary Priedmann Robertson Walker (PRW) behavior takes place throughout the early uni- 
verse. At a redshift Zcard 1 • the two terms inside the bracket become equal, and henceforth 
the second term dominates. Once the second term dominates, it drives the universe to ac- 
celerate. The energy density at which the two terms become equal is pcard = Po(l + -^Card)^, 
where po is the matter density today. The MP Cardassian model of Eq.(2) depends on three 
parameters: the numbers n and q and the density pcard- The latter can be traded for the 
observed matter mass density fij^^ (see Eq. [12] below). 

The original power law Cardassian model gave the same distance-redshift relation as 
a quintessence model with constant equation of state parameter Wg = n — 1. Generalized 
Cardassian models, on the other hand, give predictions for the distance-redshift relation 
that can be very different from generic quintessence models. Por example, some Cardassian 



^The name "Modified Polytropic" arises in the context of treating this model as a fluid; then the rela- 
tionship between energy density and pressure is roughly polytropic (see Gondolo & Preese (2002)). 
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models can satisfy the weak energy condition w > —1 even with a dark energy component 
that has an effective equation of state wx = Px/px < —1- Note that an effective wx < —1 is 
consistent with recent CMB and large scale structure data (Schuecker et al. 2002; Mclchiorri 
et al. 2002). In this paper we explore these differences and their testability for the MP 
Cardassian model. 

The Cardassian model also has the attractive feature that matter alone is sufficient to 
provide a fiat geometry. Because of the extra term on the right hand side of the Friedmann 
equation, the critical mass density necessary to have a flat universe can be modified, e.g. to 
0.3 of the usual value. Hence the matter mass density can have exactly this new critical value 
and satisfy all the observational constraints such as given by the baryon cluster fraction and 
the galaxy power spectrum. 

General Notation: 

In a fiat universe, the total energy density of the universe can be written as 

Ptotaliz) = Pcold [^^^'^(1 + Zf + Qxfxiz)] . (5) 

Here, Q'^^ is the observed matter density of the universe; we will take Q'^^ = 0.3 as our 
fiducial value. The critical density of the universe (for the ordinary Friedmann equation) 
is Pcold = 2 X 10~^^/iq gm/cm^, where ho is the Hubble constant today in units of 100 
km/s/Mpc. We take 

fix = 1 - (6) 

(for a fiat universe with total energy density iltot = 1), and fx{z = 0) = 1. The subscript 
X refers to any component of the universe that provides an additional term in Einstein's 
equation; generically it is called "dark energy" , but in the Cardassian case it is an additional 
matter term. The dark energy density is 

Px{z) = px{0)fx{z) = Pcoid^x fx{z). (7) 

If the dark energy density corresponds to a cosmological constant, then one finds that fx{z) — 
1 at all redshifts z. 

Note that for generalized Cardassian models, both terms in Eq. (5) come from matter; 

fi^* = n:^^ + fix = 1. (8) 

For the MP generalized Cardassian models of Eq. (2), the dimensionless dark energy density 
fx{z) is given by 

. . ^ Px{z)_ ^ n^^ji + zy r 
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For fixed Jl^*, it depends on the two dimensionless parameters n and q. 



In Cardassian models, the observed matter mass density fraction today is given by the 
ratio of the critical mass density of the Cardassian universe, Pc,Card — Po, ^-nd that of the 
standard universe, Pc,oid = SHq/ISttG). So the observed matter mass fraction today in the 
MP Cardassian model is (Preese & Lewis 2002) 

Pcold [l + (l + ^card)=^^(^-")]'/'' 

Inversely, we can express -Zcard; and pcard; in terms of ^2^'^ 
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Outline: First, we will compare the Modified Polytropic Cardassian model of Eq.(2) 
with existing data from supernovae and cosmic microwave background data. 

Next, we explore how plausible future SNe la data can be optimally used to constrain 
dark energy models, and whether generalized Cardassian models can be differentiated from 
generic models of quintessence and models with a cosmological constant. 

Most of this work was performed in September 2002 when all the authors were at the 
Kavli Institute for Theoretical Physics in Santa Barbara. Subsequent dispersal of all the 
authors to different parts of the country caused the conclusion of the paper to take a long 
time. 



2. Compcirison of MP Cardassian with Current Data 

In this section, we compare the Modified Polytropic Cardassian Model of Eq. (2) with 
current supernovae and cosmic microwave background (CMB) data. We will see that the 
existing data can be well fit for several choices of the parameters n and q. 

In a smooth Friedmann-Robertson- Walker (FRW) universe, the metric is given by ds^ = 
dt'^ — a^(t)[(ir^/(l — kr'^) + r'^{d6'^ + sin^ ^ o?0^)], where a(t) is the cosmic scale factor, and k 
is the global curvature parameter. The comoving distance r is given by (Weinberg 1972) 

r{z)^cH^'^, K^\n,\'/\ (13) 
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n^:^t%^x.F) = j^z'^^, (14) 

E{z') = pt^l + z'f + fx{z') + ^1.(1 + z'f] , (15) 
where = 1 - f)^" - f^x, and 

S{x) — smh(a;), Jl^ > 

— X, Qjfc = 

= sm(x), Qfe < 0. (16) 

The angular diameter distance is given by dA^z) = r{z)/ {1 + z), and the luminosity distance 
is given by dL{z) = (1 + z)r{z). 

The distance modulus for a standard candle at redshift z is 

,,,{z) ^ m - M = 5 log (^^^ + 25, (17) 

where m and M are the apparent and absolute magnitudes of the standard candle, and ^^(z) 
is its luminosity distance. 

Type la supernovae (SNe la) are our best candidates for cosmological standard candles, 
because they can be cahbrated to have small scatters in their peak luminosity (Phillips 1993; 
Riess, Press, & Kirshner 1995). 

Fig. 1 shows the measured distance modulus (actually the deviation of the distant 
modulus with respect to the expected values for an open universe with fl^^ = 0.3 and 
Qa = 0) for flux-averaged^ (Wang 2000b) SNe la data (Riess et al. 1998; Perlmutter et al. 
1999) as a function of redshift. For comparison, superposed on the data points are the 



''Here we briefly describe fliix-averaging. Due to the inhomogeneous distribution of matter in our universe, 
the hght from perfect standard candles (which ah have exactly the ik luminosity) at a given redshift 

z will experience different amounts of bending (due to matter inhomogeneity along different lines of sight) 
before reaching the observer. Hence even perfect standard candles will be observed to have a non-Gaussian 
spread in peak luminosity due to gravitational lensing (Prieman 1997; Wambsganss et al. 1997; Holz & Wald 
1998; Metcalf & Silk 1999; Wang 1999; Wang, Holz, & Mimshi 2002; Munshi & Wang 2003). If this effect 
is not properly taken into account, the estimated cosmological parameters will be biased, i.e., the estimated 
mean of the parameters will deviate from the true value of the parameters. Fortunately, the total number of 
photons from all the standard candles at redshift z should remain unchanged in the presence of gravitational 
lensing (which only redistributes the photons by bending the light from each standard candle) ; therefore the 
average peak luminosity of all the standard candles at z should be the same as the peak luminosity of a 
standard candle at z without gravitational lensing. This is the basic idea behind flux-averaging of type la 
supernova data. For details, see Wang (2000b). 
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predictions for several familiar cosmological models (dotted curves, from top to bottom): 
(n^*, Qa) = (0.3, 0.7); (0.3, 0); and (1, 0). In addition, three examples of modified polytropic 
Cardassian models from Eq. (2) are shown, all with Q'^'^ = 0.3. The three models have 
parameters n = 0.2, q = 1 (solid curve); n = 0.2, q = 2 (short-dashed curve); and n = 0.2, 
q = 3 (long-dashed curve). Note that the solid curve is equivalent to a quintessence model 
with Wq = —0.8. Also shown in Fig.l is a quintessence model with Wq{z) = — 1 + 0.5^; 
for comparison (dot-dashed hne). All three generalized Cardassian models shown satisfy 
current constraints from SNe la. So we conclude that MP Cardassian models fit the existing 
supernovae data very well. 

Fig. 2 shows the (n, q) parameter space with constraints from current observational 
data from both supernovae and the cosmic microwave background (CMB), at a fixed value 
ftrn — 0.3. The MP Cardassian model with parameters (n, q) is compared with a fiducial 
ACDM (cold dark matter) model flrn = 0.3, Qb = 0.05, Qa = 0.7, and h = 0.65. The 
constraints are derived by requiring that the MP Cardassian models agree with the fiducial 
ACDM model to within la of the measurement uncertainties of the Wilkinson Microwave 
Anisotropy Probe (WMAP) CMB data (Bennett et al. 2003) and the current SN la data 
[(Riess et al. 1998; Perlmutter et al. 1999), flux-averaged with Az = 0.05 (Wang 2000b)]. 

The region between the thick solid lines in Fig. 2 corresponds to (n, q) values of the 
MP Cardassian models that satisfy the current SN la data [fiux-averaged with Az = 0.05] 
within one sigma of the fiducial ACDM model, i.e., Ax^ = Xmpc ~ Xacdm ~ 1- Note that 
the fiducial ACDM model {Qm = 0.3, Q\ = 0.7, h = 0.65) corresponds to n = 0, q = 1. 
Complementary plots showing constraints from the power law Cardassian model with fixed 
q — 1 and varying n and fi^ are given in Sen & Sen (2002); Zhu & Fujimoto (2002, 2003). 

The dashed line in Fig. 2 indicates the constraints from WMAP CMB data. When 
the cosmological parameters are varied, the shift in the whole CMB angular spectrum is 
determined by the shift parameter (Bond, Efstathiou, & Tegmark 1997; Melchiorri et al. 
2002; Odman at el. 2002) 



where r{zdec) denotes the comoving distance to the decoupling surface in a flat universe. The 
results from WMAP data require Hor{zdec) — 3.330 ^o Jgg and hence give the allowed range 
in the shift parameter. The corresponding constraints on MP Cardassian model parameters 
(n, q) are represented by the dashed hne in Fig. 2. All models in Figure 2 below the dashed 
line lie within this allowed range for the shift parameter. We conclude that the Modified 
Polytropic Cardassian models of Eq.(2) are compatible with current supernova and CMB 
data. 




(18) 
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3. Compcirison of Models Using Simulated Future Data 

In this section, we will construct simulated type la supernova (SN la) data for three dark 
energy models: a Modified Polytropic Cardassian model, a cosmological constant model, and 
a quintessence model. We then investigate if we can recover the original theory from the 
simulated data. In particular we want to see if Cardassian cosmology can be differentiated 
from generic quintessence models or a cosmological constant by analysis of upcoming SN la 
data. We will begin with a modified polytropic Cardassian model of Eq.(2), choose specific 
values of the parameters q and n, and see how many SNe la we would expect. 

The measured distance modulus for a SN la (labeled "Z" ) is 

+ (19) 

where jj.p'' is the theoretical prediction [see Eq.(17)], and e*^') is the uncertainty in the mea- 
surement, including observational errors and intrinsic scatters in the SN la absolute mag- 
nitudes. In the simulated data set, we take the dispersion in SN la peak luminosity to be 
^iTT'int — 0.16 mag (this is the rms variance of e^'^). 

Many (but not all) models of dark energy can be characterized by an equation of state 
Wx{z) = Px{z)/ Px{z), where px{z) is the pressure. Most authors have concentrated on 
constraining the equation of state wx of the dark energy from SN data. However, it was 
shown by Maor, Brustein, & Steinhardt (2001); Barger & Marfatia (2001) that it is extremely 
hard to constrain wx using SN data. Instead, Wang & Garnavich (2001) emphasized that 
it is easier to extract information on the dark energy density px{z), instead of wx{z), from 
the data. This is because there are multiple integrals relating wx{z) to the luminosity 
distance diiz) of SN, which results in a "smearing" that obscures the difference between 
different wx{z). It is better to use px{z) directly, as it is related to the time derivative of 
the comoving distance to SN la, r'{z); hence it is less affected by the smearing effect. The 
advantage of measuring px{z) over measuring wx{z) was confirmed by Tegmark (2002). In 
their work, Wang & Garnavich (2001) assumed that p'x{z) > 0, a condition equivalent to 
the weak energy condition for those cases in which the ordinary Ftiedmann equation applies. 
Here, on the other hand, we make no such assumption. In fact, we will show that it is 
possible to determine the sign of p'xi^)- 
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3.1. Determining the sign of the time dependence of dark energy density, 

with prior on 

We simulate data for three models (see Table 1): (1) a cosmological constant model; (2) 
a MP Cardassian model with n = 0.2 and q = 2, which has p'x{z) < 0; (3) a quintessence 
model with wx{z) = — 1 + 0.5z, which has p'x{z) > 0. 



Table 1 
Dark Energy Models 



Model 


model parameters 


pW 


ACDM 

MP Cardassian model 
quintessence model 


rim = 0.3, Oa = 0.7 

= 0.3, n = 0.2, q = 2 
= 0.3, Wc^(z) = — 1 + 0.5 z 


p'xiz)<0 
P'xi^) > 



Note that MP Cardassian models can have either p'x{z) > or p'x{z) < 0. On the 
other hand, popular quintessence models in which the quintessence field tracks the matter 
field have p'x{z) > (Barger & Marfatia 2001). Fig. 3 shows the sign of p'x{z) in the (n, q) 
parameter space at fixed flm — 0.3 for MP Cardassian models [see Eq.(2)] with — 0.3. 
The arrows indicate the regions in which p'x{z) > and p'x{z) < respectively. The region 
in between indicates models with dark energy densities that are not monotonic functions of 
time. 

Note that the weak energy condition requires that the total equation of state w > —1. 
In the context of treating Cardassian models as a fiuid (Gondolo & Preese 2002), this is 

-^1 '^-") (20) 

which satisfies w > —1 for the parameter choices we are interested in: fl^^ < 1, n < 1, and 
q > 1. Therefore, all viable MP Cardassian models satisfy the weak energy condition but 
can have Wx < —1- Note an effective Wx < — 1 is consistent with recent CMB and large 
scale structure data (Schuecker et al. 2002; Melchiorri et al. 2002). 

Some scalar field dark energy models with wx < —1 have been studied previously 
(Caldwell 1999); models which are stable despite violating variants of the weak energy con- 
ditions have been found to be difficult to construct (Carroll, Hoffman, & Trodden 2003). 
Our proposal here is a different alternative to the models previously studied. 
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We will now assume that the data set is given, either from the simulated data sets 
described above, or, in the future from SNAP^. We will show that investigation of the data 
set can reproduce the sign of the time-dependence of the dark energy density, assuming one 
knows the matter density to an accuracy of 10%. 

Given our data set, we now proceed as though we have no information on where it comes 
from; i.e., we proceed as though we did not know which model it came from. Wc parametrize 
the dark energy density in order to allow us to compare it to the data set. Wc take px{z) 
to be an arbitrary function. To approximate the function, we parametrize it by its value at 
nun equally spaced rcdshift values, z-i, i = \.2,...,nun-i Zn^.^ = Zmax- The value of px{z) at 
other redshifts arc given by linear interpolation, i.e., 

Px{z) = ( — — ) Pi-i + ( ) Pi, Zi^i < Z <Zi, 

\Zi — Zi-l / \Zi — Zi-iJ 

Zq = 0, ^n^in ~ Zmax- (21) 

The values of the dark energy density pi {i = 1,2, ...,nhin) are the independent variables to 
be estimated from data; note that the number of independent variables is ribin- Again, we 
proceed as though we had absolutely no information on the function px{z), and treat it as 
a completely arbitrary function. 

The complete set of parameters, then, is 

s=(Q^^^ Pi, and nun), (22) 

where i — 1, ■■■,n}nn ^ described above. Hence our number of parameters is N — nun + 2. 
We will vary the number of bins riinn between 1 and 10, and look for the optimal fit to the 
data. To illustrate, an arbitrary function may become a good approximation to the data for 
4 bins whereas it is a miserable fit for 3 bins. 

We expand the adaptive iteration method developed in Wang & Garnavich (2001) and 
Wang & Lovelace (2001); unhke what is done in those papers, we do not restrict ourselves 
to cases where p'x{z) > 0. 

We can now determine a best fit to the set of parameters s by using a statistic, with 
(Riess et al. 1998) 



^p^p\zi\s) - p^\zi) 



2 



(s) = E^ -2 (23) 



^Note that the current SNAP design is sustantiaUy improved than before (Tarle et al. 2002). Here we 
assume that SNAP will obtain all SNe la in its survey fields up to ^ = 1.7, similar to a supernova pencil 
beam survey (Wang 2000a; Wang & Lovelace 2001). 
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where ijip\zi\s) is the prediction for the distance modulus at redshift zi, given the set of 
parameters s. Here ai is the dispersion of the measured distance modulus due to intrinsic 
and observational uncertainties in SN la peak luminosity. 

To reduce the computation time, we can integrate over the Hubble constant Hq analyt- 
ically, and define a modified statistic, with 



x' = xt--^^{Ci + - InlOj - 21n h\ (24) 

where h* is a fiducial value of the dimensionless Hubble constant h, 

At; = Hp{h = h*) = 42.384 - 5 log r + 5 log [Hor{l + z)] , (25) 



and 



I ^ ' 



a = 



It is straightforward to check that the derivative of with respect to h* is zero; hence our 
results are independent of the choice of h* . We take h* — 0.65. 

For a given choice of Ti^jn, we can minimize the modified statistic of Eq.(24) to find 
the best fit and px{z) (parametrized by pi, i = 1, 2, nbin). We can find one sigma 
error bars by finding values with (Ax)^ — 1 from the minimum. 

For each model in Table 1, we obtain four sets of best fit parameters. We apply four 
different constraints to the arbitrary function px{z) in order to discover which one allows a 
good fit. The four constraints are: 

(i) pxiz) — Px(0) = constant; i.e., a cosmological constant model; 

(ii) P'xiz) > 0; 

(iii) p'^iz) < 0; 

and (iv) completely unconstrained px{z)- 

For each of these constraints, we find the best fit parameters. 

Figure 4 shows our results: panels (a) and (b) correspond to the MP Cardassian and 
quintessence models described in Table 1 respectively. For simulated data of Model 1 of Table 
1, a cosmological constant model, we find that is estimated correctly to 1% accuracy. 
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regardless of the assumption made about p'x{z)- For Model 2 (MP Cardassian) and Model 
3 (quintessence), Fig.4(a) and (b) show the best fit under all of the four constraints 
above, for nun values ranging from 1 to 10. We find that assuming the wrong sign for p'x{z) 
leads to an estimated VL^n that differs from the assumed Vf^^ by more than 10%. The different 
constraints on the sign of p'x{z) are represented by different point types. The solid horizontal 
line is our fiducial value of VLm = 0.3 (i.e., we are assuming that that this is the true value of 
the matter density), and the dot-dashed horizontal hnes indicate 10% error bars about his 
fiducial value. We are assuming that il^ is known to within 10% from other data sets. 

These plots arc not intended to emphasize the dependence of Q"^'^ on riun- Indeed, as 
discussed above, the reason that we have found the best fit fim for a variety of Ubin values is 
simply that the parametrization of the arbitrary function px{z) may be poor for one value 
of riiyin but excellent for another; we take a given model to be a good one if it lies within the 
10% range on for several values of nun- 

Fig. 4(a) and (b) show estimated as function of nun for Model 2 (MP Cardassian 
model) and Model 3 (a quintessence model). In Fig. 4 (a), only the VL"^^ values estimated 
assuming that p'xi^) ^ consistently (i.e., for most values of nun) lie within 10% of the true 
value of = 0.3. Hence, we have indeed recovered the correct general time-dependence of 
the model underlying this set of simulated data. In Fig.4(b), only the $7°^^ values estimated 
assuming that p'x{z) > deviate by less than 10% from the true value of Vt°J^^ = 0.3. Again, 
we have recovered the correct general timc-dcpcndcncc of the model underlying this set of 
simulated data. For all three models, we have been able to correctly ascertain the sign of 
p'x{z) with this technique. 

This indicates that the estimated (for a variety of values of ribin), together with a 
10% accurate prior on 1^^'^, can be used to determine the general time-dependence of the 
dimensionless dark energy density; i.e., the sign of p'xi^)- 

3.2. Estimating dark energy density from data 

In this section we create a large number of Monte Carlo samples to see how well we can 
reconstruct the entire function px{z)- To create these samples, we need first to identify the 
best-fit model as follows. 

For a given data set, we choose the best-fit model with nun [the number of parameters 
used to parametrize the dimensionless dark energy density px{z)\ that satisfies three condi- 
tions: 

(1) it corresponds to an estimated VL^a value that deviates less than 10% from the true value; 
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(2) as we decrease from a large value, say, — 10, it minimizes the per degree 
of freedom, Xpdf = X^/i^data — t^), without significantly shifting the estimated value of flm 
(Wang & Lovelace 2001). N^ata is the number of SNe la, and u is the number of parameters 
estimated from data. Long dashed lines in Fig.4 show Xpdf function of on an arbitrary 
scale^; 

(3) if p'x{z) ^ 0, then Ubin > 1- 

Now we apply the above conditions to Fig. 4(a) and (b). In Fig. 4(a), as we decrease 
Ubin from Ubin = 10, the significant shifting in the estimated Qm occurs at Ubin = 6, which 
also has the smallest Xpdf ^'^^ ^ — ''^bin < 10. We find that for the MP Cardassian model, 
SNAP data yield an optimal Ubin = 6. In Fig.4(b), as we decrease Ubin from Ubin = 10, the 
significant shifting in the estimated flm occurs at Ubin = 3, which also has the smallest Xpdf 
for 3 < Ubin < 10. Hence for the quintessence model, SNAP data yield an optimal Ubin — 3. 

To derive the error distribution of estimated parameters ^2°^* and (i = 1, 2, Ubin, 
see Eq.(21)), we create 10^ Monte Carlo samples by adding dispersion in peak luminosity 
of Arriint — 0.16 mag to the distance modulus Hp{z) [see Eq.(17)] predicted by the best- 
fit model (i.e., assuming that the best-fit model is the true model). This is equivalent to 
making 10^ new "observations", each similar to the original data set (Press et al. 1992). The 
same analysis used to obtain the best-fit model from the data is performed on each Monte 
Carlo sample. We use the distribution of the resultant estimates of the parameters (0°^* 
and Pi) to derive the mean and 68.3% and 99.73% confidence level intervals of the estimated 
parameters. Wang & Lovelace (2001) showed that such a Monte Carlo analysis gives less 
biased estimates of parameters than a maximum likelihood analysis, i.e., the Monte Carlo 
mean of estimated parameters deviate less from the true values of the parameters. 

Fig. 5 shows the estimated dimensionless dark energy density px{z) for a generalized 
Cardassian model with n = 0.2 and q = 2 from simulated SN data from SNAP assuming 
that we know Qm to 10% accuracy. The solid line indicates the underlying true model for 
Px{z)- The horizontal dashed line near the top of the figure indicates px{z) = px{0) = 
constant, a cosmological constant model. The horizontal dotted line near the bottom of 
the figure indicates px{z) — 0. We impose px{z) > 0. The estimated fl^ values (obtained 
by reconstructing the model from the Monte Carlo samples) are listed at the bottom of 
the plot with 68.3% confidence level intervals. Where the actual value of Qm for the fake 
data set was flm = 0.3, we see that the reproduced Qm from our Monte Carlo study is 
Qjn = 0.298(— 0.023, -1-0.024). Indeed this reproduced value lies within 10% of the correct 



The scale is adjusted for the curve to fit in the figure. 
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The error bars of the reproduced estimates of px{z) have been computed using 10^ 
Monte Carlo random samples derived from the simulated data. The solid error bars and 
the dotted error bars indicate the 68.3% and 99.73% confidence level intervals respectively. 
Hence, from Fig. 5, we sec that a MP Cardassian model with a set of parameters that fit 
the current observational data, fi^* = 0.3, n = 0.2, q = 2 , can be differentiated from a 
cosmological constant model at 99.73% confidence level. We have shown the accuracy with 
which one can reconstruct the form of the dark energy density. We see that SNAP can 
indeed differentiate between different models. 



4. Discussion and Conclusion 

We have compared a particular form of Cardassian model, the Modified Polytropic 
Cardassian model of Eq. (2), with existing data from supernovae and cosmic microwave 
background measurements. We have found that current data constrain the parameter space 
of the MP Cardassian model. 

We have shown that future type la supernova (SN la) data from SNAP can differentiate 
various dark energy models (cosmological constant, quintessence, and generalized Cardassian 
expansion), assuming that is known to 10% accuracy. We have found the interesting 
result that the sign of the time dependence of the dark energy density can be determined by 
SNAP. 

Further, we have performed Monte Carlo samples to illustrate how well one can repro- 
duce the form of the dark energy density with SNAP. For example, a MP Cardassian model 
with a set of parameters that fit the current observational data, — 0.3, n — 0.2, q — 2 , 
can be differentiated from a cosmological constant model at 99.73% confidence level. 

We wish to remark on another test of generalized Cardassian models. There are two 
independent motivations for these models: 1) they may arise as a consequence of imbedding 
our observable universe as a 3-brane in higher dimensions (see, e.g., Chung & Freese (2000)), 
and 2) these models may be described in terms of a fluid interpretation, in which the dark 
energy density may be due to self-interaction of dark matter particles (Gondolo & Freese 
2002). In this second interpretation, we have developed a fully relativistic treatment of the 
resultant modified Euler's equations, Poisson equations, and energy conservation. Then we 
analyzed (Gondolo & Freese 2002) the linear growth of density fluctuations in the fluid in- 
terpretation. A more complete study of perturbation growth is in progress. Of particular 
interest is the study of the Integrated Sachs Wolfe effect in the Cosmic Microwave Back- 
ground. It is possible that the deficit of power on large angular scales (low order multipoles) 
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may be explained in generalized Cardassian models. 
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Fig. 1. — Examples of MP Cardassian models [see Eq.(2)] that satisfy current observational 

constraints from type la supcrnovac (SNe la) data. Three MP Cardassian models are shown, 
all with = 0.3: n = 0.2, q = I (solid curve); n = 0.2, q = 2 (short-dashed curve); and 
n = 0.2, q = 3 (long-dashed curve). Note that the solid curve is equivalent to a quintessence 
model with Wq = —0.8. The dot-dashed curve is a quintessence model with Wq = —1 + 0.5z. 
The dotted curves show several familiar cosmological models for comparison (from top to 
bottom): {Q^, Q^) = (0.3, 0.7); (0.3, 0); and (1, 0). 

Fig. 2. — The parameter space of (n, q) showing constraints from SNe la and CMB at fixed 
flm = 0.3. All models below the dashed line arc in agreement with the shift parameter TZ as 
measured by WMAP. All models between the solid lines are in agreement with the Type la 
SN data. 

Fig. 3. — The parameter space of (n,q) for MP Cardassian models [see Eq.(2)] with = 
0.3. The arrows indicate the regions in which p'xi^) > and p'x{z) < respectively. 
The region in between indicates models with dark energy densities that are not monotonic 
functions of time. 

Fig. 4. — This figure shows that we can indeed determine the sign of the dependence of the 
dark energy density Pxi^)] i-c, we can determine if it is increasing, decreasing, or constant 
in time. The axes show estimated fi^* as function of from the simulated data for SNAP 
for (a) a MP Cardassian model with n = 0.2 and g = 2 so that p'xi^) < and (b) a 
quintessence model with Wx{z) = — 1 -|- 0.52; so that p'x{z) > 0. The horizontal dot-dashed 
lines correspond to a 10% uncertainty on fij^*' = 0.3 ± 0.03. The different curves show the 
results obtained assuming a variety of constraints on the time dependence p'x{z) as labeled. 
The long-dashed curve shows Xpdf on an arbitrary scale. By requiring the results to lie within 
the dot-dashed lines, we recover the sign of the time-dependence of px{z)- 

Fig. 5. — Estimated dimensionless dark energy density px{z) for simulated SN data from 
SNAP, assuming that we know VLm to ~10% accuracy [and the correct general time de- 
pendence of the dark energy density, see discussion in Section 3. A]. The error bars of the 
estimated px{z) have been computed from 10^ Monte Carlo random samples derived from 
the simulated data. The sofid error bars and the dotted error bars indicate the 68.3% and 
99.73% confidence level intervals respectively. The reproduced estimates of values are 
fisted at the bottom of the plot with 68.3% confidence level intervals. 



-20- 



1.5 



1 - 



0.5 - 







-0.5 - 



-1 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

MP Cardassian models quintessence models 
n°^=.3. n=.2 ^^(z) = -1+0.5Z 

q=l 

q=2 

q=3 





Hi 



(.3..7) 

(■3,0) 



_ flux-averaged SN la data (Wang 2000b), 
_ 92 SNe la (Riess et al. 1999 & 

Perlmutter et al.l999) 
J I I I I I I I I I I i_ 



(1.0) 



J L 







0.5 



1 

z 



1.5 



Fig. 1. — Examples of MP Cardassian models [see Eq.(2)] that satisfy current observational 
constraints from type la supernovae (SNe la) data. Three MP Cardassian models are shown, 
all with = 0.3: n = 0.2, q = 1 (solid curve); n = 0.2, q = 2 (short-dashed curve); and 
n — 0.2, q = 3 (long-dashed curve). Note that the solid curve is equivalent to a quintessence 
model with Wq — —0.8. The dot-dashed curve is a quintessence model with Wq — —1 + 0.5z. 
The dotted curves show several famihar cosmological models for comparison (from top to 
bottom): (Q^, Qa) = (0.3, 0.7); (0.3, 0); and (1, 0). 
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Fig. 2. — The parameter space of (n, q) showing constraints from SNe la and CMB at fixed 
Qm = 0.3. All models below the dashed line are in agreement with the shift parameter TZ as 
measured by WMAP. All models between the solid lines are in agreement with the Type la 
SN data. 
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Fig. 3. — The parameter space of (n,q) for MP Cardassian models [see Eq.(2)] with = 
0.3. The arrows indicate the regions in which p'xi^) > and p'xi^) < respectively. 
The region in between indicates models with dark energy densities that are not monotonic 
functions of time. 
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Fig. 4. — This figure shows that we can indeed determine the sign of the dependence of the 
dark energy density Pxi^)] i-G-, we can determine if it is increasing, decreasing, or constant 
in time. The axes show estimated Q^* as function of riinn from the simulated data for SNAP 
for (a) a MP Cardassian model with n — 0.2 and g = 2 so that p'x{z) < and (b) a 
quintessence model with wx{z) = — 1 + 0.5z so that p'x{z) > 0. The horizontal dot-dashed 
lines correspond to a 10% uncertainty on fl'^^ = 0.3 ± 0.03. The different curves show the 
results obtained assuming a variety of constraints on the time dependence p'x{z) as labeled. 
The long-dashed curve shows Xpdf '^^ arbitrary scale. By requiring the results to lie within 
the dot-dashed lines, we recover the sign of the time-dependence of px{z). 



(2300 SNe la with 4m|„,=.16) 
- MP Cardassian model. 
0^=0.3, n=.g, q=2 




■ ; 99.7% C.L. 

- 

10* Monte Carlo samples 
n„ = 0.298 (0.275, 0.322) 



0.5 1 1.5 2 

Z 



-24- 



Fig. 5. — Estimated dimensionless dark energy density px{z) for simulated SN data from 
SNAP, assuming that wc know to ~10% accuracy [and the correct general time de- 
pendence of the dark energy density, see discussion in Section 3. A]. The error bars of the 
estimated pxi^) have been computed from 10^ Monte Carlo random samples derived from 
the simulated data. The solid error bars and the dotted error bars indicate the 68.3% and 
99.73% confidence level intervals respectively. The reproduced estimates of flm values are 
hsted at the bottom of the plot with 68.3% confidence level intervals. 



